Space-time dispersion of graphene conductivity 
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We present an analytic calculation of the conductivity of pure graphene as a function of frequency 
io, wave-vector k, and temperature for the range where the energies related to all these parameters 
are small in comparison with the band parameter 7 = 3 eV. The simple asymptotic expressions 
are given in various limiting cases. For instance, the conductivity for kvo <C T <C uo is equal to 
a(u), k) = e 2 /Ah and independent of the band structure parameters 7 and vo- Our results are also 
used to explain the known dependence of the graphite conductivity on temperature and pressure. 
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INTRODUCTION 

Recently the properties of graphene attracted special 
attention (see, for instance |lj and references therein). 
The matter of fact the wide variety of carbon materials 
such as graphite, graphene tubules, and fullerenes con- 
sists of graphene (i.e., a single layer of graphite) sheets 
shaped in different manner. 

Two-dimensional graphene has a very simple band 
structure, which can be obtained with help of the symme- 
try consideration or in tight-binding approximation Q. 
It was shown in Refs. [3, 0] that the energy bands of 
graphene are degenerated at the corners of the 2D Bril- 
louin zone K = (0,47r/3V3a), where a=1.44 A is the 
interatomic distance (see Fig. 1). This is the Dirac-type 
spectrum but it is massless and two-dimensional. It was 
demonstrated that due to the symmetry arguments such 
gapless spectrum with the conic point in the 3D case 
turns out stable with respect to the Coulomb interaction 
0. One can proof that this stability remains also for 
the 2D graphene spectrum with the conic point. Such 
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FIG. 1: The Brillouin zone of graphene and the electron spec- 
trum in the vicinity of the K points 

simple band structure can be used in analytic calcula- 
tions of various thermodynamic and transport properties 
of graphene. An example of such calculations was pre- 
sented in Ref. where the imaginary part of the di- 
electric function Ime(w) of graphite has been calculated. 
The gapless band structure of graphene also results in 
the unusual behavior of conductivity. Despite the enor- 



mous number of paper devoted to carbon materials, the 
analytical expression of graphene conductivity was not 
yet derived to our best knowledge. 

The standard diagrammatic approach to calculation 
of transport properties of impure metals is restricted by 
the fundamental requirement that the mean free path of 
carriers t = v$t must be much larger than the electron 
wavelength A = h/pp, i.e. ipp 3> 1. This condition 
cannot be evidently satisfied in the case of graphene when 
the Fermi surface degenerates to the points. We avoid 
this difficulty addressing the problem when the graphene 
sample is clean enough but temperature is finite (TV 3> 
I). In this case the temperature appears instead of Fermi 
energy and electrons obey Boltzmann statistics. 

In the present paper, we calculate analytically the fre- 
quency and wave- vector dependence of the graphene con- 
ductivity <Jij(u>,k) at finite temperature T for the rela- 
tively low frequencies u> <C 3 eV, when only the nearest 
7r bands can be taken into consideration. 



HAMILTONIAN AND CURRENT 

Let us remind briefly the properties of the graphene 
electron spectrum. The tight-binding approximation al- 
lows to write down the following effective Hamiltonian 
2x2 matrix 

where h(p) = 7 {e ip * a + 2e~ ip " a / 2 cos ( Py aV3/2)} and 
7 =< ip(a, 0)\Ho\ip(0, 0) > is the ppir transfer integral, 
which is of the order of 3 eV. The dispersion law of 
graphene, ex, 2 (p) can be obtained by means of its di- 
agonalization: 

£1,2 (p) = ±7 jl + 4cos (3p x a/2) cos (p y aV3/2)+ 

4cos 2 (p y aV3/2)| 1/2 . (2) 
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In the vicinity of the K point, the matrix elements of the 
effective Hamiltonian can be expanded as 

MP) = v o{iPx +Py), 

where vq — i^a/2. This results in the linear dispersion 
relation £1.2 (p) = ±«oP- Calculating the conductivity, we 
can use this expansion for to <C 7 « 3 eV. Let us pass to 
the calculation of the electric current, following the paper 
by Abrikosov [t| , where the dielectric function of Bi-type 
metals was considered. We assume that the external 
field is described by the vector potential Aj . The current 
operator has the form 



(3) 

where x' — > x, v x , x and (m 1 ) x i x are the velocity and 
mass operators correspondingly, the tilde in the notation 
ip means that the operator is taken in the interaction 
representation with 



V = -- c \ ip+ix^vi^^Aiix^x. 



(4) 



We will calculate the current in the linear approxima- 
tion in Ai(x). Therefore, the second term in Eq. © can 
be taken in zeroth approximation with respect to the in- 
teraction Expanding the first term in Eq. ((3J) to the 
first order in A4 (x) , we get the retarded correlator of four 
ip operators which has to be averaged over the Gibbs en- 
semble. At finite temperatures, the Fourier component 
(with respect to the spatial coordinates and imaginary 
time) of this correlator 
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(5) 



is expressed in terms of the temperature Green's func- 
tions: 



G(p) 



H(p)Y 



In Eq. JSJ the notations p±are used, the summa- 
tion is carried out over the fermionic frequencies u> n = 
2irT(n + l/2), while the trace operation is performed 
over the band index of the Hamiltonian H(p). The lat- 
ter can be easily carried out in the representation where 
the Hamiltonian has diagonal form: 

Tr {v l Qv j G} = VixGnviiGu + ^22^22^22^22 



functions one finds 



T^2Gii{p + )G 2 2(p-) 
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e 2 (p+) +ei(p_) 
1 

£i(P-) iw+-£ 2 (p+) 

= /oMp-)]-/oMp + )] 

iui - £ 2 (p+) + ei(P-) ' 



(6) 



where we have taken into account that the photon fre- 
quencies u>i = 2-nlT are "even" in the Matsubara tech- 
nique; /o is the Fermi distribution function. Analytical 
continuation of the expressions similar to Eq. JJJ) into the 
upper half plane of the complex frequency can be per- 
formed by simple substitution iloi —> lu + i5 with 5 — ► 0. 

Let us notice that the current has to vanish when the 
vector potential does not vary in time. Since the second 
term in Eq. @ is time independent, one can omit it, 
subtracting from the first term its value at u> = 0. In 
result 



<Tij(w, k) 



E 



ie 



(7) 



d>V{/o[e„(p-)]-/o[e a (p + )]} 



„ [eo(p+)-Eo(P-)][w-£o(p+)+Eo(P-)] 
dVi2^2i{/o[ei(P-)] ~ /oMp+)]} 



-2uj 



Mp+) - ei(p_)] V - [£ 2 (P+) - ei(P-)] 2 } 



This expression acquired the factor 4 due to summation 
over spin and over six points of the K type (two per each 
Brillouin zone). 

Hitherto we did not use any peculiarities of the 
graphene spectrum besides the number of bands. Thus 
Eq. J7J has a general character. For graphene, the matrix 
of the velocity near the point K in the band representa- 
tion is determined by the Hamiltonian Q: 



Wo 











(8) 



where are unit vectors along the coordinate axis direc- 
tions. The unitary transformation which transforms the 
Hamiltonian from the band representation to the diago- 
nal one has the form 



U 



1 f Py~ iPx Py - ipx 
V2p \ ~P P 



Its application to the velocity matrix gives: 



+vl 2 G22V 3 21 Gii + v l 21 Gnv{ 2 G22- 

Then one can perform the summation over u> n in stan- 
dard way. For instance, for the cross product of the Green 



jj-i v jj _ ^0 ( e x p x + e y p y i(e x p y - e y p x ) 



P \ i(^xPy QyPx) ^xPx ^yPy 



The first term in Eq. J7J corresponds to the intra-band 
electron-photon scattering processes. In the limit of the 
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high carriers concentration kvo <C (T, Ep), it leads to the 
usual Drude expression 



intra / , r\\ 
xx ( w , °) = 



El a p 



(iLU — T 1 ) ^— 
y 1 a=l 



d 2 p 2 df [e a (p)) 



de 



(9) 

with 1/r — > (t is the transport scattering time). The 
second term in Eq. J7J owes its origin to the inter-band 
electron transitions. The real part of this contribution 
(let us recall, that iuji — > lu + iS) at k — > is reduced to 
the well-known expression for the absorbed energy due 
to the direct inter-band transitions. 



In the limit of low temperatures T << lu, the energies 
e ~ T « oj are essential and we can take Taylor of the 
integrand over T /lu : 



9T 

- — [In2 + 6C(3)(1» 2 ], 

LU Z 



(12) 



where £(3) = 1.20. One can see that the term with In 2 
exactly cancels the intra-band contribution l]lUp. 

In the limit of high temperatures T >> oj, the lead- 
ing contribution to the integral (| 1 1 11 originates from the 
region lu < e < T. One can obtain with the logarithmic 
accuracy: 



ASYMPTOTIC BEHAVIOR OF CONDUCTIVITY 



Ip = - 



1 f T de 



8T 



- = -— ln- 

e 8T n uo' 



(13) 



Let us pass to the discussion of the the pure graphcnc 
conductivity in absence of gate voltage, when the chem- 
ical potential is equal to zero. The integral J7J) can be 
performed analytically for various limiting cases. 



Small spatial dispersion kvo Cu,T 



Collecting the Eqs. l[n jl .lfT2" ]) . and JEJ, one can write: 



al n * er {u, k) = — tanh — - — x (14) 

(T/cj)[ln2 + 6C(3)(T/w) 2 ], kv « T « uj, 
(w/16T) In (T/lj), kv « u « T. 



Putting k = and integrating over angle, one can 
find that the off-diagonal elements of conductivity vanish, 
while the diagonal ones are equal to: 



a xx (cj,0) = <7 yy {u},0) 
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ede 



dfo(e) 



ds 



fo(-e) - f (e) 



de / jo 
The first (intra-band) contribution 

ie 2 T 



4e 2 



l (oJ,k) = 2 In 2 



kvn « lu, T 



(10) 



(we write explicitly the Planck constant in the final ex- 
pressions). This result was obtained in Ref. Q using the 
Drude expression, Eq. (EJ. 

The inter-band term in Eq. 1)10(1 has the form: 



"(^,0) = — tanh — 



4T 



(11) 



ie 2 cu 



P 

7T J 



de e 
^4? tanh 2T' 



where P denotes the operation of taking the principal 
value of an integral. This integral can be presented in 
more convenient form 



P 



ds 



io u 2 - 4e 
-2P 



tanh 



e 

2T 

de 



1 



o 



(w 2 -4e 2 )(exp (e/T) + l)' 



Large spatial dispersion lu <C kvo,T 

Now let us consider the limit of large dispersion 
kvo >> lu. We choose the direction of k along the x- 
axis. In the intra-band term of Eq. J7J, one can expand 



£2,i(P±) = \J e 2 (p) + (vok/2) 2 ± pkv 2 cos <p « (15) 
£ (p)±pkvl cos 4>/2e , 



where £o (p) — \J £ 2 (p) + (vok/2) 2 , e (p) = v p, and (f> is 
the angle between k and p. 

Now we evaluate the integral in the Eq. 10 for 
cr ra™( w ! First of all, let us note that changing of 
the variable of integration p — > p results in 

[LU - £ a (p+) + £ a (p_)] _1 + [LU + £ a (P+) - £a(P-)] _1 

2lu 



LU 



[e a (p+) - £a(p-)] 2 ' 



Using the expansion (|15|l one can see that in the signif- 
icant domain of integration, determined from the condi- 
tion lu <C pkvQ /eo (p) , the integrand of intra-band part 
cr™* ra in Eq. (0, can be presented as follows 



> dj o(£o) 



Luen 



x deo {pkvo cos < 



where v x — vo cos 4>. Thus, cos <f> disappears in the in- 
tegrand. But for the intra-band contribution in the 
transversal & y n v tra component of conductivity, we have 
v y = v o sin </>, and therefore, the longitudinal and 
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transversal components (respect with the k - direction) 
are different: 



intra ( , , i \ 
xx ( w , fe ) = 



{2T/kv ) 2 \n2, 



ijj < kv < T, 



In (2-\/ kvoT I to) exp (—kvo/2T), lu < T < kvo, 



(4 T/fc« )l n2, uj<kv <T, 
\J 7rfcwo/Tcxp (—kvo/2T), u><T<kvQ. 



What concerns the inter-band contribution, here we 
can put £2,1 = ±e (p)- Taking into account that vf 2 — 
ivQ sin cj) and integrating over <fr, one can obtain 



'(v,k) 



-in 



deo Eq 
-2-tanh— . 

kvo/2 £ a Z1 



Evaluating this integral, we obtain for the inter-band 
contribution 



(i/4r 

1/kvo, 



(1/4T) ln(4T/fcu ), lj < kv < T, 
uj < T < kvg. 



CONCLUSIONS 



The expression Eq. (|14|l allows us to estimate the con- 
ductivity of pure graphene, t _1 = 0. In this case, the 
first term plays the leading role (for small k and at low 
temperatures, when tanhw/T = 1: 



<t(uj, k) = — , kv Q <T<w. 



(16) 



Let us underline that this conductivity results from the 
electron transitions between two intersecting bands at 
the K points of the Brillouin zone. Remarkable fact is 
that its value turns out to be universal, independent of 
any parameter of graphene, like 70 or vq. 

Very recently the electric field effect in graphene was 
investigated experimentally Q]. Indeed, the universal 
conductivity behavior for the samples with different car- 
rier concentrations was found. In spite of fact that the 
experiment was performed in conditions to <C T, different 
from Eq. I|16[l . the minimal value of measured conduc- 
tivity cr( min ) = 2e 2 /nh was found close to our prediction. 

At high temperatures, when the condition kvo «w« 
T is fulfilled, the conductivity Eq. 1)11) ll becomes imagi- 
nary: 



a(u,k) = 2iln2 



e 2 T 
ttTiuj 



and it depends on temperature. 

The behavior of such type can be observed on experi- 
ments involving the plasmon modes. The dispersion law 
of such modes for 2D systems is gapless, 



to 



hi k . 



2e 2 7Tn2 

Jj2 2 ' 



where £oo is the lattice dielectric constant. 

So far we considered the 2D graphene sheet. The re- 
sults obtained can be immediately applied to the 3D 
graphite if one neglect the interaction between the layers. 
Then the integration with respect to the p z component 
of the quasi-momentum in the Brillouin zone gives just 
the additional factor l/c 2 in comparison with the con- 
ductivity of graphene, where c z is the distance between 
the layers in the z-direction. For instance, Eqs. JHJ and 
(|10[) acquire the factor l/c z . In the low- frequency limit 
lu -C 1/t, conductivity can be estimated as 



r (3) 



(u,k -> 0) = 21n2 



e 2 Tr 
nhcz 



The number of phonons in the 2D graphene at low tem- 
peratures (T <C To , To is the Debye temperature) is pro- 
portional to T 2 . Since the Fermi surface is assumed to be 
small (ef < T < To), all scattering angles are essential 
in this region of temperatures. For the electron collision 
rate, which is determined by the electron-phonon inter- 
action, one can write t _1 = aT 2 /Tjj with the constant a 
of the order of unity. Thus the in-layer resistivity turns 
out to be linear in temperature: 

nhac z T aT ^ , 

^2^1^- 32 °^^ Cm - W 

According to Ref. &\, p = 60 /xficm at T = 300 K 
in agreement with the above estimate (To ~ 2000-ftT for 
graphite, a ~ 1). The equation l|17fl answers the question 
discussed in Refs. 0,0 on the pressure dependence of the 
graphite resistivity. We see that the resistivity decreases 
under the pressure because the inter-layer distance c z 
decreases and the Debye temperature To grows. 

This work was supported by the RFBR, Grant No 04- 
02-17087. The authors are grateful to D.V. Livanov for 
valuable discussions. 



[1] K.S. Novoselov et al., Nature, 438,197 (2005). 
[2] R. Saito, G. Dresselhaus, and M.S. Dresselhaus, Physical 
Properties of Carbon Nanotubes (Imperial College Press, 
London, 1998). 
[3] P.R. Wallace, Phys. Rev. 71, 622 (1947). 
[4] J.C. Slonczewski and P.R. Weiss, Phys. Rev 99, 636(A) 
(1955). 

A. A. Abrikosov, S.D. Beneslavsky, Sov. Phys. JETP 32, 
699 (1971). 



[5 



■5 



[6] T.G. Pedersen, Phys. Rev. B 67, 113106 (2003). 
[7] A. A. Abrikosov, Sov. Phys. JETP 17, 1396 (1963). 
[8] L. Edman, B. Sundqvist, E. McRae, E. Litvin-Staszewska, 
Phys. Rev. B 57, 6227 (1998). 



[9] R.G. Arkhipov, V.V. Kechin, A.I. Likhter, Yu. A. 
Pospelov, Sov. Phys. JETP 17 1321 (1963). 



